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Abstract —In this work we present new distributed controllers 
for secondary frequency and voltage control in islanded micro- 
grids. Inspired by techniques from cooperative control, the pro¬ 
posed controllers use localized information and nearest-neighbor 
communication to collectively perform secondary control actions. 
The frequency controller rapidly regulates the microgrid fre¬ 
quency to its nominal value while maintaining active power 
sharing among the distributed generators. l\ming of the voltage 
controller provides a simple and intuitive trade-off between the 
conflicting goals of voltage regulation and reactive power sharing. 
Our designs require no knowledge of the microgrid topology, 
impedances or loads. The distributed architecture allows for 
flexibility and redundancy, and eliminates the need for a central 
microgrid controller. We provide a voltage stability analysis 
and present extensive experimental results validating our designs, 
verifying robust performance under communication failure and 
during plug-and-play operation. 

Index Terms —Microgrid, distributed control, secondary con¬ 
trol, inverters, voltage control 


1. Introduction 

E conomic factors, environmental concerns, and the rapidly 
expanding integration of small-scale renewable energy 
sources are all pushing the incumbent centralized power gener¬ 
ation paradigm towards a more distributed future. As a bridge 
between high-voltage transmission and low-voltage distributed 
generation (DG), the concept of a micro grid continues to gain 
popularity [l]-[5]. Microgrids are low-voltage electrical dis¬ 
tribution networks, heterogeneously composed of distributed 
generation, storage, load, and managed autonomously from 
the larger primary network. Microgrids can connect to a 
larger power system through a Point of Common Coupling 
(PCC), but are also able to “island” themselves and operate 
independently [2]. Islanded operation of a microgrid could be 
planned, or could occur spontaneously if a fault triggers the 
disconnection of the microgrid from the primary grid. 
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Fig. 1: Low-detail schematic of microgrid control architecture. 


Energy generation within a microgrid can be quite heteroge¬ 
neous, including photovoltaic, wind, micro-turbines, etc. Such 
sources generate either DC power or variable frequency AC 
power, and are interfaced with a synchronous AC grid via 
power electronic inverters. It is through these inverters that 
cooperative actions must be taken to ensure synchronization, 
voltage regulation, power balance and load sharing in the 
network [6]. Control strategies ranging from centralized to 
completely decentralized have been proposed to address these 
challenges [7], and have subsequently been aggregated into a 
hierarchical control architecture [2] (Figure 1). 

The control hierarchy consists of three levels. The first 
and most basic level is primary control, which stabilizes the 
microgrid and establishes power sharing. Although centralized 
architectures have been used for primary control [7], in order 
to enhance redundancy and enable plug-and-play functionality, 
the current standard is to employ proportional control loops 
locally at each inverter. While successful for stabilization, 
these decentralized “droop” controllers force the bus voltages 
and the steady-state network frequency to deviate from their 
nominal values [6], [8], [9]. 

This leads naturally to the next level in the hierarchy, termed 
secondary control. Broadly speaking, the goal of secondary 
control is to remove the aforementioned deviations in both 
global frequency and local voltage [6]. Centralized techniques 
for secondary control have been well studied in high-voltage 
transmission and distribution networks [10]. These centralized 
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strategies have also been applied in the context of microgrids, 
and the term “secondary” has been broadened to include 
additional control goals such as reactive power sharing [11], 
[12], harmonic compensation, and voltage unbalance [2], [6], 
[8], [13]. The final level of tertiary control is concerned with 
global economic dispatch over the network, and depends on 
current energy markets and prices. 

Several recent works (see Section II-C) have proposed 
secondary control strategies for microgrids. However, to date 
no single control strategy has been able to offer a flexible, 
plug-and-play architecture guaranteeing frequency and voltage 
regulation while maintaining precise active and reactive power 
sharing among non-identical DGs. Currently, this combination 
of goals appears infeasible with decentralized control using 
only local information (voltage, power, ect.) at each DC [14]. 
As such, communication between DGs has been identified as 
a key ingredient in achieving these goals while avoiding a 
centralized control architecture [12], [14]-[19]. 

In this paper we build on our previous theoretical and 
experimental works [15], [16] and introduce a general and 
fully distributed framework for secondary frequency and volt¬ 
age control in islanded microgrids. Our designs overcome the 
limitations of existing strategies by combining decentralized 
proportional droop control and integral control with distributed 
averaging algorithms. We therefore refer to our proposed con¬ 
trollers as DAPI (Distributed Averaging Proportional Integral) 
controllers. These controllers use decentralized control actions 
and sparse communication among neighboring DC units to 
achieve precise frequency regulation, active power sharing, 
and a tunable trade-off between voltage regulation and reactive 
power sharing. The distributed architecture eliminates the need 
for a central supervisory control: additional DGs are integrated 
through a low-bandwidth communication link to an existing 
DG, with the communication topology being a tunable design 
variable. The DAPI controllers are model-free, in the sense that 
they require no a priori knowledge of the microgrid topology, 
line impedances or load demands. 

There are four main technical contributions in this paper. 
First, In Section III we highlight and clearly demonstrate 
a fundamental limitation of voltage control: precise voltage 
regulation and precise reactive power sharing are conflicting 
objectives. The presentation frames and motivates our sub¬ 
sequent controller designs. Second, in Section IV we review 
the frequency DAPI controller [15] and introduce the voltage 
DAPI controller. This new voltage DAPI controller accounts 
for the conflict between voltage regulation and reactive power 
sharing by allowing for a tunable compromise between the two 
objectives. We build intuition for our design by detailing sev¬ 
eral tuning strategies. Taken together, the two DAPI controllers 
form a distributed duo for plug-and-play microgrid control. 
Third, in Section V we present a small-signal voltage stability 
analysis of the microgrid under DAPI control, derive sufficient 
conditions on the controller gains and microgrid parameters for 
closed-loop stability, and study the transient performance of 
the system under changes in the controller gains. Finally, in 
Section VI we present extensive experimental results validat¬ 
ing our DAPI designs. The experimental microgrid consists 
of four heterogeneous DGs in a non-parallel configuration. 


with high R/X connections and distributed load. We validate 
our designs, and move beyond our theoretical results by 
demonstrating controller performance under communication 
link failures and plug-and-play operation. 

Section II contains a review of standard control strategies 
for microgrids. In particular. Section II-B reviews primary 
droop control, while Section II-C provides a detailed review 
of secondary control for both frequency and voltage. Our main 
contributions are housed in Sections III-VI, with concluding 
remarks being offered in Section VII. 

II. Review of Microgrid Control 
A. Problem Setup and Review of Power Flow 

In this work we consider microgrids consisting of n buses 
which are either DGs or loads. For inductive lines of reactance 
Xij connecting bus i to bus j, the active and reactive power 
injections Pi and Qi at bus i are given by [10] 

A = y'. (la) 

Xij 

Q* = 5 - E” 1 ^ ’ (lb) 

Xi ^J=1 Xij 

where Ei (resp. Oi) is the voltage magnitude (resp. voltage 
phase angle) at bus i and Xi = l/(X]j=i 


B. Review of Primary Droop Control 

A complete survey of primary control is beyond the scope 
of this paper; we provide here a short summary. The objective 
of primary (droop) control is to stabilize the network and 
establish a proportional sharing of load among the DGs. In 
islanded operation, inverters are controlled as grid-forming 
Voltage Source Inverters, having controlled frequencies and 
voltage magnitudes [3]. To accomplish this, a foundation of 
control loops must be established to regulate the current, 
voltage, and output impedance of the inverter. This is achieved 
via an inner control loop for the current, an outer control loop 
for the voltage, and a virtual impedance loop ensuring the 
output impedance is desirable at the line frequency (lower half 
of Figure 7) [3], [20]. The reference inputs for the voltage 
loop are supplied by droop controllers, which are heuristic 
controllers based on active/reactive power decoupling [2], [6], 
[8], [18], [21]. For inductive lines, the controllers specify the 
inverter frequencies uji and voltage magnitudes Ei by^ 


- rriiPi, 

(2a) 


(2b) 


where cc* (resp. F^*) is a nominal network frequency (resp. 
voltage), and Pi (resp. Qf) is the measured active (resp. 
reactive) power injection. The gains rrii^ni are the droop co¬ 
efficients. In [15] a large-signal stability analysis of (la)-(2a) 
was completed, yielding the steady-state network frequency 


<^ss — cc* + 


Po 

E n 5 

i=l rrii 


(3) 


^Without loss of generality, in islanded mode we consider the droop 
equations (2) without power set points; if desired these can be included as 
iUi = iU TTli(^Pi -Pz,set) Snd Ei = E (^Lset)- 
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where Pq is the total active power load in the microgrid. 
Note that the steady-state frequency (3) is different from 
the nominal cc*. Large-signal stability analysis of the voltage 
droop controller ( 2 b) remains an open problem; see [ 22 ] for 
results on a related droop controller. For non-inductive lines, 
the appropriate droop controllers take other forms [3]. 

C. Review of Secondary Control 

The removal of the steady-state frequency and voltage 
deviations generated by the droop controllers ( 2 a)-( 2 b) is 
accomplished by “secondary” integral controllers. 

1) Frequency Regulation: Many techniques have been sug¬ 
gested to restore the network frequency, ranging on the spec¬ 
trum from centralized to decentralized [ 6 ], and each with its 
own advantages and disadvantages. One centralized technique 
is to mimic Automatic Generation Control from bulk power 
systems. This is implemented using area control errors on 
slow time-scales, a centralized integral controller, and one- 
to-all communication [10]. However, this centralized approach 
conflicts with the microgrid paradigm of distributed generation 
and autonomous management. A decentralized technique is to 
use a slower integral control locally at each inverter [ 8 ]. This 
implicitly assumes that the measured local frequency is equal 
to the steady-state network frequency, and therefore relies on 
a separation of time-scales between the fast, synchronization- 
enforcing primary droop controller and the slower, secondary 
integral controller [ 8 ], [9]. Except in special cases, this ap¬ 
proach destroys the power sharing properties established by 
primary control [14], and is too slow to dynamically regulate 
the grid frequency during rapid load changes. 

In [12], [23] control strategies were proposed in which 
DG units communicate their frequencies, voltages and reactive 
power injections to one another in order to perform secondary 
control and share active and reactive power. The methods have 
two drawbacks: first, all inverters must communicate with all 
other inverters, requiring a dense communication architecture. 
Second, the controller gains must be finely tuned in order to 
maintain active power sharing; see [19] for a detailed analysis. 

2) Voltage Regulation: In high-voltage networks, the shar¬ 
ing among generators of reactive power demand is usually not 
a major concern due to capacitive compensation of both loads 
and transmission lines; voltages at generators are therefore 
regulated to fixed values by the excitation system [10]. Voltage 
regulation has subsequently been adopted as the standard for 
voltage secondary control in microgrids [2], [ 8 ]. However, 
in small-scale microgrid applications, the low ratings of DG 
units, small electrical distances between units, and the lack of 
static compensation requires an accurate sharing of reactive 
power demand among DGs to prevent overloading. In Section 
III we highlight how voltage regulation and reactive power 
sharing are conflicting objectives. 

Due to the line impedance effect, the voltage droop con¬ 
troller ( 2 b) is unable to share reactive power demand among 
even identical inverters operating in parallel [ 6 ]. In [18], an 
alternative primary droop controller was proposed for reactive 
power sharing among parallel inverters with the same rated 



Fig. 2: Schematic of DGs operating in a parallel microgrid. 


voltages. The method requires each unit to have a measure¬ 
ment of the common load voltage, limiting its applicability 
in complex microgrid scenarios. Similarly, the centralized 
secondary control architecture proposed in [ 11 ] for reactive 
power sharing requires each unit to communicate with a 
central controller. The distributed voltage controller proposed 
in [12], [23] require all DGs to communicate with all others 
directly. Moreover, since the controller regulates DG voltages 
to their nominal values, it is be unable to share reactive power 
between heterogeneous units connected through varying line 
impedances. See [11], [12] and the references therein for more. 

III. Fundamental Limitations of Voltage Control 

In this section we illustrate the fundamental conflict between 
two secondary control goals: voltage regulation and reactive 
power sharing. For simplicity of exposition, we focus on a 
parallel microgrid consisting of two identical DGs connected 
to a common distribution bus (Figure 2). The reactances of the 
two lines connecting the DGs to the common bus are different; 
in particular, Xqi > X 02 . 

Figure 3 depicts the E-Q droop characteristics before and 
after a standard, voltage-regulating secondary control action. 
Without secondary control, the inverters operate at voltages Ei 
and E 2 with reactive power injections Qi and Q 2 (solid black 
line). Since Qi 7 ^ Q 2 , reactive power is not shared; this is 
the “line impedance effect”. Application of voltage-regulating 
secondary control ensures that both DG voltage magnitudes 
are restored to the common rating E* (dotted blue and green 
lines are the post-secondary control droop characteristics). 
Note however that the inverter power injections change to 
Qi < Qi and Q 2 > Q 2 - The application of standard 
secondary control therefore worsens the already poor sharing 
of reactive power between the DGs. 

For the same problem setup. Figure 4 depicts the E-Q 
droop characteristics before and after a power sharing en¬ 
forcing secondary control action is taken.^ While the identical 
inverters now proportionally share the reactive power by both 
injecting Q", the resulting voltage values E'{ and E 2 are more 
dissimilar than they were with only primary control. 

Table I collects the relationships between voltage magni¬ 
tudes and reactive power injections for the different control 

^This control action is not uniquely determined; there are many shiftings 
of the droop characteristics which lead to power sharing (Section IV-C). 
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^ Conventional Droop (DG 1 & 2) 



Fig. 3: E-Q droop and standard secondary control for two parallel inverters 
with identical ratings, operating through reactive lines with Xqi > Xo 2 - 



Fig. 4: E-Q droop and power sharing secondary control for two parallel 
inverters with identical ratings, operating through reactive lines with 
Xoi > Xo2. 


TABLE I: RELATIONSHIPS BETWEEN VOLTAGE MAGNITUDES AND 
REACTIVE POWER INJECTIONS FOR DIFFERENT CONTROL 


ACTIONS 

Control Method 

Voltage Magnitudes 

Reactive Powers 

Primary Control 
Standard Sec. (') 
Power Sharing {") 

E2 < El < E* 
E'^=E[= E^ 

No Relationship 

Qi ^ Q2 

Q'l < Qi < Q2 < Q2 
Q’i = Q’i = Q” 


actions described above. We observe that — except under 
special circumstances — precise voltage regulation leads to 
large errors in reactive power sharing, as shown in Figure 3. 
Conversely, the objective of reactive power sharing does not 
uniquely determine the DG bus voltages, and when imple¬ 
mented improperly can result in poor voltage profiles as shown 
in Figure 4. The accuracy of reactive power sharing that can be 
achieved therefore depends on both the upper and lower limits 
for the DG voltage magnitudes, and on the homogeneity of the 
line reactances. We conclude that an ideal secondary voltage 
controller should allow for a tunable compromise between 
voltage regulation and reactive power sharing. 

IV. Distributed Averaging Proportional Integral 
(DAPI) Controllers for Microgrids 

As mentioned in Section I, communication has been identi¬ 
fied as an essential ingredient for high-performance secondary 


control. We now introduce the DAPI controllers, which com¬ 
bine droop and integral control with averaging algorithms from 
multi-agent systems [24]. To build intuition for our designs, we 
first briefly review continuous-time averaging over networks. 


A. Review of Continuous-Time Distributed Averaging 

The communication layer between DG’s will be described 
by a weighted graph G(V, 5, A) where V = {1, ..., n} is a 
labeling of the DGs, f C V x V is the set of communication 
links, and A is the n x n weighted adjacency matrix of 
the graph, with elements = Oji > 0. In particular, one 
writes that (i,j) G £ if node i sends information directly 
to node j, and in this case Oji > 0. Thus, the sparsity 
pattern of the adjacency matrix A encodes the topology of the 
communication layer (Figure 5). If to each node i G {1, ..., n} 



( ^ 

^12 

0 

aiA 

CLI2 

0 

Ci23 

024 

0 

Ci23 

0 

Os4 

\ai4 

024 

a34 

0 ) 


Fig. 5: Example of adjacency matrix construction for four DGs, with 
V = {1, 2, 3, 4} and S: = {(1, 2), (2,1), (2, 3), (3, 2), (3, 4),...}. 


we assign a scalar value Xi, a commonly studied update rule 
is for node i to adjust its value Xi according to 

E tl 

Xj). (4) 

Equation (4) is called continuous-time distributed averaging, 
or “consensus”. To interpret (4), define the convex weights^ 
Wij = CLij/{Y^]^=i ^ik)^ and rearrange (4) to obtain 

1 . 

Z^y=l ^ij 

Thus, with time-constant 1/ variable Xi evolves 

toward a weighted average of its neighbors values Xj, 
with averaging weights wij. If the communication network 
G(V,f,A) is connected, this dynamic process results in all 
variables Xi converging to a common value Xi = Xj = const, 
(see Remark 1) [24], [25]. We now apply these ideas from 
continuous-time distributed averaging to microgrid control. 


B. Frequency Regulation and Active Power Sharing 

We propose the distributed-averaging proportional-integral 
(DAPI) frequency controller 


UJi — CC* — TTliPi -\- Cli , 

dl2 ^ 

hi ^ ^ ^ ) 

i=i 


( 6 a) 

( 6 b) 


where Vti is the secondary control variable and hi is a positive 
gain. The first equation ( 6 a) is the standard droop controller 
with the additional secondary control input ^2^. To understand 
the second equation, we consider two cases. 


^The weights are convex because Wij > 0 and Wij = 1. 
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Conventional Droop (DG 1 & 2) 



Fig. 6: Droop characteristics with (dashed lines) and without (solid lines) 
secondary control action. One can interpret the DAPI secondary control 
action as a uniform shifting of all droop characteristics by an amount 

(jJ CUss • 


Case 1 (A = 0): In this case, there is no communication 
among DGs, and (7^ is the integral of the local frequency 
error — cc* with gain 1/ki. In steady-state, the derivative 
on the left of (6b) is zero, and hence uji = uj* for each DG i. 
That is, the network frequency has been regulated. However, 
depending on the initial conditions and controller gains, the 
variables Qi may converge to different values, and shift their 
respective droop curves by different amounts. This unwanted 
degree of freedom leads to poor active power sharing. 

Case 2 (A 7 ^ 0): Now, consider the case where we include 
the diffusive averaging terms aij{fti — ftj). As before, in 
steady-state the derivative on the left-hand side of (6b) must 
be zero, and hence uji = cc*. However, from the discussion in 
Section IV-A, we also must have Qi = Qj for all inverters i, j. 
That is, the DGs must agree on how much to shift the droop 
characteristics. This ensures that all droop curves are shifted 
by the same amount equal to cc* — ccgs (Figure 6), guaranteeing 
active power sharing is maintained. This performance is not 
dependent on the controller gains ki and , which determine 
only the transient behavior of the controller (see Table II). 

Remark 1: (Communication Requirements for DAPI 
Control). DAPI control requires that neighboring DG units ex¬ 
change information to collectively perform secondary control. 
To ensure power sharing among all units, the communication 
network among DGs must be connected: there must be a path 
in the communication graph between any two nodes, as in 
Figure 5. While here we consider the controller in continuous¬ 
time with bi-directional communication, our assumptions can 
be relaxed to allow for asymmetric, asynchronous and discrete¬ 
time communication with delays [24], [25]. 


C. Voltage Regulation and Reactive Power Sharing 

As noted in Section II-C, the E-Q droop controller (2b) 
is unable to share reactive power between DGs. Moreover, 
in Section III we described the conflict between voltage 
regulation and reactive power sharing. With these problems 


in mind, we propose the second DAPI controller 


Ei = — riiQi + Ci , 

i=i 




(7a) 

(7b) 


where is the secondary control variable, (5* is the ith DGs 
reactive power rating, and /Si, Hi are positive gains. The n x n 
matrix B with elements hij > 0 is the adjacency matrix of 
a communication network between the DGs. The secondary 
controller (7b) achieves a tunable compromise between voltage 
regulation and reactive power sharing. We consider four cases: 

Case 1 (P = 0, B ^ 0): In this case the first term in 
(7b) is disappears, leaving only the second term. Steady- 
state requires the derivative on the left-hand side of (7b) 
to be zero, which occurs if and only if QijQX = Q^jQ] 
for all inverters. Thus, the steady-state is a power sharing 
configuration. The secondary control variables converge to 
values which shift the individual droop curves as necessary to 
establish proportional power sharing, see Figure 4. However, 
as discussed in Section III, under such a control action DG 
voltages can drift quite far from their nominal values. 

Case 2 (/3 7 ^ 0, S = 0): In this case the second term in 
(7b) disappears, and the controller reduces to the standard 
decentralized voltage-regulating secondary control discussed 
in Section III. Reactive power is shared poorly (Figure 3). 

Case 3 (/3 7 ^ 0, S 7 ^ 0): In this regime (7a)-(7b) achieves 
a compromise between reactive power sharing and voltage 
regulation based on the relative sizes of the gains (Si and hij. 

Case 4 (Smart Tuning): As a specialization of Case 3, 
consider having a specific DG i implement the controller (7b) 
with Pi 0 and hij = 0, while the all other DGs j ^ i 
implement (7b) with Pj = 0 and hjk / 0."^ That is, DG i 
regulates its voltage to the nominal value, and the voltages at 
DGs j i are then controlled to share power in a manner 
consistent with the voltage regulation of DG i (cf. Section 
III). This tuning sets up a “leader-follower” [26] relationship 
among the DGs, where the voltages at DGs j ^ i will form a 
cluster around the voltage value of Ei = E* of DG i. 

The above cases are tested experimentally in Section VI-A. 

Remark 2: (Remarks on DAPI Control). The communica¬ 
tion layers between DG units described the adjacency matrices 
A and B are design variables of the DAPI controllers. This 
customizable architecture allows for design flexibility. For ex¬ 
ample, to add redundancy against communication channels be¬ 
ing permanently disconnected, supplementary communication 
can be introduced. Note that the communication architecture 
need not mirror the electrical topology of the network (Figure 
9), and that the controllers do not rely on high-gain techniques 
such as feedback linearization [27]. A detailed schematic of 
the DAPI control architecture is shown in Figure 7. 

The time-constants ki and Hi in (6b) and (7b) allow for a 
precise tuning of the secondary control speed. A conventional 
choice is to make ki and Hi sufficiently large, enforcing a time- 
scale separation between primary and secondary control. This 


"^This directed communication tuning requires that DG i sends information 
to at least one neighbor. 
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Low Bandwidth 



partial analysis for a simpler controller can be found in [28]. 
In Section V-A we present a small-signal stability analysis 
of (7a)-(7b), along with sufficient conditions which ensure 
stable operation. In Section V-B we explore the effect of the 
controller gains in (6),(7) on the transient performance of the 
closed-loop system. 

A. Small-Signal Stability of Voltage DAPI Control 

To avoid unnecessary technical complications, we model 
any delay in adjusting the output voltage in (7a) with a 
simple low-pass filter, yielding the dynamic system dEi/dt = 
— {Ei — E^) — riiQi -f Ci, and assume loads are impedances 
collocated with DGs. Both of these assumptions can be relaxed 
at the expense of more complicated formulae. Under the 
standard decoupling assumption in which reactive power is 
related strongly to differences in voltage magnitudes [10], the 
reactive power injection (lb) at the ith DG takes the form 

Qi = + Ei ~ Ej) , (8) 

where Ijoad is diagonal matrix of load susceptances and 
^bus = i^ microgrid’s bus admittance matrix [10]. 
In vector notation, the system equations (7),(8) are^ 

E = -{E -E^)-NQEe, (9a) 

Ke = -^{E-E*)-L^[Q*]-^Q, (9b) 

Q = [E]YE, (9c) 


Fig. 7: Block diagram of proposed control architecture for a single DG. For 
simplicity we have have abbreviated di = Yfj=i ^ij - 


TABLE II: QUALITATIVE EFFECTS OF CONTROLLER GAINS 


Gain | Qualitative Change Upon Increase 


ki 

Ki 

aij 

bij 

di 


Slows frequency regulation at DG i 

Slows voltage regulation / Q-sharing at DG i 

Faster P-sharing between DGs i and j 

Improved steady-state Q-sharing between DGs i and j 

Improved steady-state voltage regulation at DG i 


however is not required - our experimental results suggest 
that primary and secondary control can be performed on 
similar time scales without stability issues or performance 
degradation.^ Table II provides a simple qualitative reference 
for the effects of the control parameters in (6) and (7). 

V. Stability & Performance of DAPI Control 

A large-signal nonlinear stability analysis of the frequency 
DAPI controller (6) can be found in [15]. While the secondary 
frequency controller (6b) will never destabilize the primary 
controller (6a), the secondary voltage controller (7b) can 
potentially destabilize (7a). This possibility exists due to the 
previously discussed conflict between reactive power sharing 
and voltage regulation. A full nonlinear stability analysis 
of the voltage/reactive power DAPI controller (7a)-(7b) is 
extremely challenging and beyond the scope of this article; a 

^In fact, in [15] it was shown that the frequency controller (6a)-(6b) is 
stabilizing for any choice of gains ki. 


where and e are the vectors of voltage magnitudes, 

voltage set points, and secondary control variables, /3 and 
K are diagonal matrices of controller gains, Q and Q* are the 
vectors of DG reactive power injections and reactive power 
ratings, Lc = diag(5^J^i bij) — B is the Laplacian matrix [25] 
corresponding to the communication network among the DGs, 
and Y = — (Ibus + ^load)- When implementing the controller 
(7) in practice, the voltages Ei will remain near their nominal 
values ^*. We can exploit this to obtain a linear dynamic 
system by making the approximation that [E] [T^*] in (9c); 

details on this approximation technique can be found in [29]. 
After making this approximation and inserting (9c) into (9b), 
the nonlinear system (9a)-(9b) becomes the linear system 

X = Wx + u , (10) 

where x = {E, e), u = (E’*, and 

/-Wi lA f -{I^+N[E*]Y) I„\ 

V-^2 On] + L,[QTHE*]Y) On] ’ 

where In (resp. On) is the n x n identity matrix (resp. zero 
matrix). For future use, we note that all eigenvalues of —Wi 
are real and negative since Wi is similar to a symmetric M- 
matrix, as can be verified by using the similarity transform 
TWiT-^ where T = Ny/‘^[E*j]-'^/'^. We now derive 
sufficient conditions under which the linearized system (10) 
is exponentially stable. Specifically, we will assume that the 
characteristic polynomial dei{sl 2 n — W) = 0 of (10) has a root 
in the closed right-half complex plane, and derive conditions 

^Here [z] denotes the diagonal matrix with the vector 2 ; along the diagonal. 
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under which this assumption is contradicted. These conditions 
will therefore ensure all characteristic roots are in the left- 
half complex plane, and thus ensure stability. Since —Wi has 
negative eigenvalues, it follows that det{sin + Wi) ^ 0, and 
using the Schur complement determinant formulae for block 
matrices we may simplify the characteristic polynomial as 

det(5/2n -W) = det{s‘^In + sWi + W 2 ) = 0 . (11) 

Since the determinant is zero, the matrix s‘^In + sWi -]- W 2 
must be singular, and therefore the polynomial (11) has a 
solution if and only if x^{s‘^In + sWi -f 1 ^ 2 )^ = 0 for 
some real vector x of unit length. The latter is simply a scalar 
quadratic equation of the form 5^ + (Ti 5 -h (T2 =0, where 
ai = x^Wix and a 2 = x^W 2 X. If it is true that 

-^min {Wi+Wl)>0, (12a) 

\^UW2+W^)>{), (12b) 

where Aniin(') is the smallest eigenvalue of the matrix argu¬ 
ment, then 0 ( 1 , 0(2 >0 and all solutions of 5 ^ + 0(i5 + 0(2 =0 
satisfy Re(5) < 0 by the Routh-Hurwitz criterion. This con¬ 
tradicts our assumption that the characteristic polynomial has 
a closed right-half plane root, and hence under the conditions 
(12) the linearized system is exponentially stable. 

Let us now physically interpret the stability conditions (12). 
The first condition (12a) restricts the DGs from being too 
dissimilar. For example, if all DGs have the same droop gains 
rii and voltage set points then Wi is a scalar values times 
Y and (12a) is always true. For dissimilar DGs, the intuition 
for (12a) is that given equal voltage set points, DGs with 
high ratings should be connected to the microgrid through stiff 
lines of high admittance. To understand the second condition 
(12b), we first consider the case of pure voltage regulation 
(Case 2 in Section IV-C) where A / 0 Lc = On- Then 
W 2 = > Ois diagonal and (12b) is satisfied. The voltage 

regulation gains pi always act to stabilize the system. Since 
eigenvalues are continuous functions of matrix parameters, 
the system is also stable for non-zero but sufficiently small 
power sharing gains hij. In the more general case where the 
power sharing gains hij are also non-zero, the condition (12b) 
properly accounts for the complicated interplay between the 
microgrid’s electrical stiffness matrix Y and the averaging 
control action Lc in the product [F^*]F. Intuitively, 

(12b) will be satisfied when all line impedances are sufficiently 
uniform and all DGs are sufficiently similar, since in this case 
reactive power sharing is not in strong conflict with the line 
impedance effect (cf. Section III). The stability conditions (12) 
are both satisfied in all experiments presented in Section VI. 

B. Transient Performance of DAPI Control 

The impact of the controller gains ki , Oij ^^i^ pi and hij on 
the steady-state equilibrium was discussed in detail in Section 
IV and summarized in Table II. We now examine the impact of 
these controller gains on the system’s transient performance. 
To do this, we consider a case study with four DGs (Figure 9), 
with the system parameters of Table III. The communication 
network among the DGs is a ring, and the controller gains 
Oij and hij are given by (13); Oij are constants, while hij 


are parameterized by a single constant h. For all DGs i = 
1,..., 4, the other control parameters ki^ Ki and Pi are taken 
as uniform constants /c, n and P, respectively. The nominal 
values for these gains are k = 1.7 s, k, = Is, P = 1.2, and 
h = 180 V (the same as in Study lc of Section VI). 

We increment the gains independently in intervals around 
their nominal values, and for each iteration we (i) numerically 
determine the system operating point from (1),(6),(7), (ii) lin¬ 
earize the closed-loop system around the operating point, and 
(iii) plot the eigenvalues of the linearization. These eigenvalue 
traces are displayed in Figure 8, where black crosses denote 
the eigenvalues for the nominal controller gains and arrows 
indicate the direction of increasing gain. 

Eigenvalues on the real axis are strongly associated with 
the frequency dynamics (6a)-(6b), while complex conjugate 
eigenvalues are associated with the voltage dynamics (7a)- 
(7b). These conjugate eigenvalues lead to an underdamped 
voltage response: physically, this is a manifestation of the 
line impedance effect, which the controller (7a)-(7b) must 
overcome to establish reactive power sharing. As the frequency 
time-constant k is increased (blue), real eigenvalues move 
towards the origin leading to slow, smooth frequency/active 
power response. Conversely, decreasing k leads to fast (but 
still overdamped) frequency regulation. Increasing Oij has 
an effect nearly identical to decreasing k, and we have 
therefore omitted the trace and held Oij constant. Increasing 
the voltage time-constant k (red) causes the underdamped 
conjugate eigenvalues to collapse onto the real axis, leading to 
an overdamped voltage/reactive power response for sufficiently 
slow secondary control. Increasing either feedback gain h or 
P (green and gold) results in an increasingly underdamped 
voltage/reactive power response. 

Taken together. Table II, the stability conditions (12), and 
the eigenvalue traces of Figure 8 provide a solid foundation for 
understanding the DAPI controllers (6)-(7). Our experimental 
results demonstrate that despite the simplifying assumptions 
in the preceding analysis, the DAPI controllers (6)-(7) can be 
tuned for both stability and high performance. 

VI. Experimental Results 

Experiments were performed at the Intelligent Microgrid 
Laboratory (Aalborg University, Denmark) to validate the 
DAPI controllers presented in Section IV. A schematic of the 
experimental setup is shown in Figure 9, consisting of four 
DGs interconnected through impedances. Loads are present 
locally at units 1 and 4, and units 1 and 4 are rated for 
twice as much power as units 1 and 3 (Table III). The DAPI 
controllers (6a)-(7b) were implemented in Simulink®, with 
measurements recorded via a dSPACE® 1006. See [2], [3] 
for details on the inner voltage, current and impedance control 
loops. 

This section is organized into four studies, beginning with a 
characterization of controller performance, and then examining 
robustness under communication link failure, heterogeneous 
controller gains, and plug-and-play operation. The commu¬ 
nication structure is shown in Figure 9, with the adjacency 
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Fig. 8: Eigenvalue traces of closed-loop system (I),(6),(7) as controller 
gains are varied. Arrows indicate the direction increasing gain. The absence 
of a trace indicates that the parameter under consideration has negligible 
effect on the respective eigenvalue. System parameters are taken from Study 
Ic of Section VI. Black crosses indicate eigenvalue locations for the nominal 
gains used in Study Ic. Several fast eigenvalues are omitted for clarity. 



Fig. 9: Schematic of the experimental microgrid setup, consisting of four 
DGs interconnected through heterogeneous impedances. Loads are 
collocated at DGs one and four. Red dotted lines denote communication 
links. 


matrices A = [a^y] in (6b) and B = [bij] in (7b) being 
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(13) 


where b > 0 varies depending on the study under consider¬ 
ation. All other parameters are as reported in Table III. All 
plots are color-coded in correspondence with Figure 9: DG 1 
(blue), DG 2 (red), DG 3 (green), and DG 4 (brown). 


A. Study 1: Controller Performance 

Studies in this section illustrate the performance of the DAPI 
controllers (6a)-(7b) under various controller tunings. In all 
four sub-studies la-Id, only primary droop controllers (6a) 
and (7a) are running up to f = 7 s, at which time the secondary 


TABLE III: ELECTRICAL AND CONTROL PARAMETERS 


Parameter 

Symbol 

Value 


Electrical Setup 

Nominal Erequency 


50 Hz 


DC Voltage 

Fdc 

650 V 


Nominal Voltages 


325.3 V (230 V rms) 

Eilter Capacitance 

C 

25 fiF 


Eilter Inductance 

Lf 

1.8 mH 


Output Impedance 

Lo 

1.8 mH 


Line Impedance (1,2) 

Zl 2 

R = 0.8n, L = 

3.6 mH 

Line Impedance (2,3) 

Z23 

R = 0An, L = 

1.8 mH 

Line Impedance (3,4) 

Z34 

R = 0.7Q, L = 

1.9 mH 


Control Parameters 

Parameter 

Symbol 

DGs I&4 

DGs 2&3 


Rated Active Power 

p* 

1400 W 

700 W 


Rated Reactive Power 

q) 

800 VAr 

400 VAr 


P — uo Droop Coeff. 

Q — E Droop Coeff. 

mi 

rii 

2.5 • 10-3 iscL 5 . 10-3 

1.5-10-3^ 3-10-3 

rad 

VAr 

Int. Erequency Gain 

ki 

I.7s 

I.7s 

Int. Voltage Gain 

Ki 

Is 

Is 



controllers (6b) and (7b) are activated. The local load at DG 
unit 4 is detached at t = 22 s, then reattached at t = 36 s. 

First considering the frequency dynamics in the top por¬ 
tion of Figure 10, the frequency deviation experienced under 
primary droop control is quickly eliminated by the DAPI 
controller (6a)-(6b), and frequency regulation is maintained 
throughout load changes with minimal transients. Active 
power is accurately shared amongs the heterogeneous DGs 
throughout the entire runtime. This robust frequency and active 
power behavior is identical in all other sub-studies, and thus 
we omit the plots due to space considerations. The remainder 
of studies la and lb in Figures 10 and 11 highlight the 
conclusions drawn in Section III regarding the limitations of 
voltage secondary control. 

a) Tuning for pure reactive power sharing: Figure 10 
shows results for the voltage DAPI controller (7a)-(7b) tuned 
for power sharing (/3i = 0, b = 50V), with no attempt at 
voltage regulation, as in Case 1 of Section IV-C. While reactive 
power is shared accurately, voltage magnitudes deviate from 
their nominal values = 325.3 V (cf. Figure 4). 

b) Tuning for pure voltage regulation: Conversely, Fig¬ 
ure 11 reports results for the same controller tuned to regulate 
voltage levels without reactive power sharing = 2.2, 
b = 0 V), as in Case 2 of Section IV-C. While voltage levels 
are tightly regulated, reactive power sharing among the units 
is poor (cf. Figure 3). As explained in Section III, the poor 
performance in Figure 11 is a general property of all voltage 
controllers strategies that exactly regulate DG voltage levels. 

c) Compromised tuning: Figure 12 displays the results 
for Study Ic, in which the DAPI controllers (6a)-(7b) are 
implemented with b = 180 V and uniform controller gains 
Pi = 1.2, as in Case 3 of Section IV-C. Considering the voltage 
dynamics, the voltage DAPI controller (7a)-(7b) achieves a 
compromise between voltage regulation and reactive power 
sharing. Voltage magnitudes are roughly clustered around 

= 325.3 V, while reactive power is approximately shared. 

d) Smart tuning: Figure 13 displays the results for Study 
Id, in which the DAPI controllers (6a)-(7b) are implemented 
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Voltage Frequency 



Active Power Injections 



DAPI Outputs (U) 



Voltage Magnitudes Reactive Power Injections DAPI Outputs (e) 




Fig. 10: Study la - reactive power sharing without voltage regulation, with control parameters 6 = 50 V, /3i = /32 = /^s = = 0. In correspondence with 

Figure 4 of Section III, the quality of voltage regulation is quite poor. 


Voltage Magnitudes Reactive Power Injections 



DAPI Outputs (e) 



Fig. 11: Study lb - voltage regulation without reactive power sharing, with parameters 6 = OV, /3i = ^2 = ^3 = ^4 = 2.2. In correspondence with Figure 
3 of Section III, the quality of reactive power sharing is quite poor. 


Voltage Magnitudes 



Reactive Power Injections DAPI Outputs (e) 




Fig. 12: Study Ic - a compromise between voltage regulation and reactive power sharing, with control parameters = [32 = [^3 = = 1-2, h = 180 V. 


with b = 100 V and ^2 = Pi = Ps = Pa = 0, in accordance 
with the discussion of Case 4 in Section IV-C. In comparison 
with the voltage dynamics of Study Ic, the voltage regulation 
in Figure 13 shows a slight improvement, while the reactive 
power sharing is noticeably improved, maintaining accurate 


sharing through load changes and during transients. Note 
that this performance improvement has been achieved while 
reducing the controller gain b which enforces reactive power 
sharing. Due to this reduction in controller gain, the ringing 
in the reactive power signal during transients is noticeably 
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Voltage Magnitudes Reactive Power Injections 




DAPI Outputs (e) 



Fig. 13: Study Id - accurate reactive power sharing and good voltage regulation, with control parameters = /34 = 0, ^2 = 4 and b = 100 V. 


improved from Study Ic to Study Id, in agreement with the 
stability and root locus analyses of Section V. 

B. Study 2: Communication Link Failure 

In this study the communication link (Figure 9) between DG 
units 3 and 4 fails at t = 2 s. At t = 7 s the local load at unit 4 
is detached, and is reattached at f = 18 s. Control parameters 
are the same as in Study Id. As the results in Figure 14 show, 
the DAPI controllers (6a)-(7b) maintain the high performance 
from Study Id despite the absence of the communication link 
between DG units 3 and 4 (cf. Remark 2). 

C Study 3: Non-Uniform Controller Gains 

We examine the behavior of the frequency DAPI controller 
( 6 a)-( 6 b) under inhomogeneous controller gains. Control pa¬ 
rameters are the same as in Study Id, except for variations 
in the integral gains ki = 1.5s, /c 2 = Is, ks = 2 s and 
k^ = 0.5 s. The results are displayed in Figure 15. Note that 
the inhomogeneous controller gains leads to varying transient 
responses for the DGs, but the steady-state behavior and 
stability of the system is unchanged. This illustrates the utility 
of the gains ki and Ki in tuning the transient response of the 
DAPI-controlled microgrid. 


Voltage Frequency Voltage Magnitudes 








■ 0 5 10 15 20 25 

Time (s) 


15 20 25 30 

Time (s) 


Fig. 14: Study 2 - DAPI performance under communication link failure. 
Control parameters are the same as in Study Id. 


D. Study 4: Plug-and-Play Functionality 

The plug-and-play functionality of the controllers was tested 
by disconnecting unit 3 at t = 10 s, and reconnecting it at 
t = 30 s. A synchronization process was used in the downtime 
to synchronize unit 3 with the remaining microgrid before 
reconnection. Control parameters are the same as in Study 
Id, and the results are displayed in Figure 16. As in previous 
experiments, the DAPI controllers (6a)-(7b) maintain accurate 
power sharing and frequency and voltage regulation before, 
during, and after the plug-and-play procedure, with minimal 
transients. The bus voltages and bus frequencies remain well 
regulated despite the disconnection of DG 3. 

VII. Conclusions 

We have introduced a general distributed control method¬ 
ology for primary/secondary control in islanded microgrids. 
By leveraging distributed averaging algorithms from multi¬ 
agent systems, the DAPI controllers achieve frequency reg¬ 
ulation while sharing active power proportionally, and can 


Voltage Frequency 



Active Power Injections 



Time (s) 


Fig. 15: Study 3 - DAPI performance with heterogeneous controller gains. 
Control parameters are the same as in Study Id. 


be tuned to achieve either voltage regulation, reactive power 
sharing, or a compromise between the two. A small-signal 
stability analysis has been presented for the voltage DAPI 
controller along with a performance study, and the controllers 
have been validated through extensive experimental testing. 

Large-signal stability of the voltage controller (7a)-(7b) 
remains an open analysis problem. Moreover, the secondary 



































































































IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS 


II 


Voltage Frequency 



Voltage Magnitudes 



Active Power Injections 



Reactive Power Injections 



Time (s) 


Fig. 16: Study 4 - DAPI performance under plug-and-pIay operation. 
Control parameters are the same as in Study Id. 


control goals of voltage regulation and reactive power sharing 
ignore an important factor for microgrid stability: the voltage 
levels at loads which are not collocated with DGs. An open 
problem is to design a controller guaranteeing that voltage 
levels at non-collocated load buses remain within tolerances 
while maintaining a high level of performance. 
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